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Abstract 

A fully microscopic theory of electron spin relaxation by the D'yakonov-Perel' 

type spin-orbit coupling is developed for a semiconductor quantum well in an 

ambient magnetic field applied perpendicular to the plane of the well. We 

derive Bloch equations for an electron spin in the well and determine explicit 

microscopic expressions for the spin relaxation times. The dependencies of 

the electron spin relaxation rates on magnetic field strength, temperature, 

and the lowest subband energy of the quantum well are analyzed. 
PACS numbers: 72.25.Rb, 73.63.Hs, 72.10.Di 
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I. INTRODUCTION 



Spin phenomenology in semiconductor structures has been at the focal point of research 
interest over the past few years in connection with proposals for spin-based quantum devices 
[1]. Spin implementation of quantum computation, optical switches, magnetic memory cells, 
etc., calls for a precise knowledge of spin dynamics and, in particular, the spin relaxation 
rates. The basic mechanisms responsible for spin relaxation are those of D'yakonov-Perel' 
(DP) [2,3], EUiott-Yafet [4], and Bir-Aronov-Pikus [5]. It was shown [6] that for III-V 
and II- VI compounds, which are the most promising materials for device purposes, the DP 
mechanism dominates at moderate temperatures and low hole concentrations. The spin 
relaxation time due to the DP mechanism has generally been expressed in the following 
semiphenomenological form [2,3]: 

for bulk semiconductors, and 

for quantum well structures, where a describes conduction band spin splitting due to lack 
of inversion symmetry [a = 0.07 for GaAs), Eg is the band gap, and T is the Kelvin temper- 
ature {kB = 1). The numerical coefficient Q depends on the orbital scattering mechanism, 
and {pD is the average square of momentum in the quantum well growth direction. It 
should be noted that these formulas involve the average momentum relaxation time, r^, as 
a phenomenological parameter. 

In the present paper we develop a fully microscopic theory of spin dynamics and ap- 
ply it to a quantum well structure in the presence of an external magnetic field, directed 
along the well growth direction, taking account of various scattering mechanisms. Our the- 
ory facilitates the microscopic determination of Tp and its temperature and magnetic field 
dependencies in terms of the fundamental material parameters. 



2 



Our analysis involves explicit recognition of the two stages in the relaxation process 
corresponding to the relaxation time hierarchies involved in (a) electron thermalization due 
to interaction with phonons, and (b) spin relaxation. In the first stage of solution, (a), 
we determine the relaxation rates and fluctuation characteristics of electron orbital motion 
due to coupling to the phonon bath. Spin relaxation dynamics (the slowest process in the 
system) can be neglected in this stage. 

The second stage, (b) , proceeds with analysis of the spin relaxation process due to spin- 
orbit interaction, wherein the orbital degrees of freedom are considered as an effective heat 
bath, having the characteristics determined in the first stage. A standard set of Bloch 
equations with two distinct relaxation times (longitudinal relaxation time, Ti, responsible 
for spin magnetic moment relaxation, and transverse relaxation time, T2, responsible for 
decoherence) is derived in this second stage. In both stages of our analyses we employ the 
method proposed in Ref. [7] and developed in Ref. [8,9]. 



II. FORMULATION 

The Hamiltonian of part (a) is given by 

H = Ho + U{r^,t) + Hpf„ (3) 
where spin and its interactions are neglected in the first stage, and 

"222 ^ ^ 

is the Hamiltonian of a two-dimensional electron in a quantum well with harmonic confine- 
ment in 2;-direction having frequency cuq. We assume the temperature to be low enough, 
so that only the lowest energy subband of the quantum well is occupied. In this case the 
motion of an electron can be described by means of two electron in-plane velocity compo- 
nents, which in the presence of a constant, uniform magnetic field directed along the z-axis 
B = (0, 0, B) are given by 
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14 = - - , Vy = -(py- -Ay) , [V^ , Vy] _ = 

m \ c / m \ c / 



ihjjr 



m 



where A = (A^, Ay, 0) is the vector potential, B = ■§^Ay — -^-Ar^, 



UJr. 



(5) 



|e| B /mc is the 



cyclotron frequency, and [...,...]_ denotes a commutator. The phonon Hamiltonian has the 
form 



(6) 



where is the phonon energy and 6^ and 6k are creation and annihilation operators, 
respectively. The term [/(rx,t) in Eq.(3) describes the electron coupling to phonons and 
impurities. 

The electron-phonon interaction is given by 



(7) 



where 



(8) 



is the phonon heat bath variable (C is the strength of the electron-phonon coupling), and 
the electron variable conjugate to this operator is defined as 



L3/2 



ikxrx(t) 



(9) 



Here, is the volume of the crystal, = {x,y), and /(/c^) = exp ^— 4^^;^^ is the electron 
confinement form factor. The response function of uncoupled phonon heat bath variables, 
fkit'iti), and their correlation function, Mk(t;ii), are given by 



(10) 



and 



i(o), 



(11) 
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where Q^\t) are the unperturbed phonon variables, Q(t — ti) is the Heavyside unit step 
function, and [...,...]_|_ denotes an anticommutator. 
The electron-impurity interaction is described by 

Ue-^{r^) = E f^k/(fc.)e^'^-'--, (12) 

k 

where are the spatial Fourier components of the impurity potential with correlation 
function 

$k=(|^[t/k,t/-kU^. (13) 

In summary, 

C/(r, t) = C/e-pft(rx, t) + Ue-i{vA_) = - J] {Q,.{t) + U^) X^{t). (14) 

k 

III. ORBITAL DYNAMICS 

Employing the Hamiltonian of Eq.(3), we obtain Langevin-like operator equations of 
motion determining in-plane electron dynamics, as given by (the derivation is presented in 
Appendix A): 

^V,{t) + uj.Vyit) + [VM. Vy{t)] = Ut). 

jVy{t) - u,V,{t) + Gy [VM, Vy{t)] = e,(t), (15) 

with fluctuation sources Cx,yit) and collision terms Gx [Vx{t)',Vy{t)] given by Eqs.(A8) and 
(A9). It should be emphasized that the expressions for fluctuation sources (Eq.(A8)) are 
obtained from microscopic analysis, and it is possible to calculate their correlation functions 
of any order. In particular, the correlation function of the fluctuation sources (in the case 
of weak coupling or Gaussian statistics of the unperturbed phonon variables) is given by 



X ( (Mk(t,ti) + <fk) ( ^ [Xk(t),X_k(ti)]+\ + i?k(r) [Xk(t),X_k(ti)]_ ) ) , (16) 



where Rk{t — ti) — h {(pk{t — ti) + (pk{ti — t)) /2i . 

We can separate the electron velocity operator into its average and fluctuation parts, 

V.,y(t)^{K,y{t))+V,,y{t). (17) 

Due to relaxation processes only the fluctuating part is nonzero, (V^,j/(i)) = 0, and the 
equations of motion for the fluctuating components take the forms 

jV^{t) + u;,Vy{t) + [VM,Vy{t)] - {G, [V,{t); Vyit)]) = Ut) 

jVyit) - uMt) + Gy [VM, Vy{t)] - {Gy [V^, Vyit)]) = ^ " ^ 

To further simphfy Eqs. (18), we have to calculate the (anti)commutators involved in 
Eqs. (16), (A8), and (A9). This procedure is presented in Appendix B resulting in the 
following simplified equations: 



{jt + ^'^^^ ~ ^"^^ ~ ^^^^^ = ^^^^ 

where the damping rate, — — ly, and the frequency shift, 8 = 5x = 5y, are given by 
+,,M„.(M.)}exp{-5^} (20) 



and 



^-^y = :^Z2 Wik,) drr | (M^ir) + $0 ^ sin (^^rj 



+^.(r)cos^^.J|cxp|-^^|. (21) 

Here, t — t — ti and Tp^(kj_) is given by Eq.(Bll). It follows from Eq. (19) that the Fourier 
transforms of velocity correlation functions are given by 

(^[V;(c^);V^/5]+^= j dre'-^ (^[Vo^{t + T),V0{t)]}l , a,(3 = x,y, (22) 



{\ [V.{oj)- V^]\ = l\ Vy]\ = ^ (- + \, , (23) 

\2 / \2 ^/ 2 V(a;-u;c) + 7o (a; + u;c) +7o/ 

and 

1 [l-.M;y,] A = - A = ^ (- - , (24) 

2 ^/ \2 ^/ 2z V(u;-a;^) +7^ {u + Uc) +7o/ 

where we neglected terms of the of order 5/uJc « 1. The correlation function of the 
fluctuation forces is defined as K{uj) = Kaaif^)-, where 

Kap{oj)= j dtTe'^^K^f,{t + r,t) 

_«,Meo.(|i.)}exp{-^}. (25) 

To carry out the r— and k— integrations in Eqs.(20) and (25), we have to employ exphcit 
expressions for the response and correlation functions for a particular scattering mechanism. 
In the present paper we consider polar optical phonons, deformational acoustic phonons and 
charged impurities as possible scattering mechanisms. Their corresponding response and 
correlation functions are listed below. 

Optical phonons: 

^k'^(^) = -p^«in(^oT)r;(r), Mk^^(r) = --^ cos (I^qt) coth(^), (26) 

where VLq is the optical phonon frequency, 1/e* = l/e^o — l/^o; (^oo and eo are the hf and 
static permittivities of the crystal, respectively). 
Deformational acoustic phonons: 

<{r) = ^ sin {ukr) rjir), M^ir) = cos (ukr) coth(^), (27) 

where D is the deformation potential , p is the crystal density, and u is the sound velocity. 
Static charged impurities have only correlation function, given by 



where ro is the screening radius, — ^aZ^, ria is the impurity concentration for species 
a, and Za is their charge number. 

We neglect the cross correlations of scattering processes and, consequently, they are 
additive in both the damping rate and in the correlation function of the fluctuation sources: 



o,OP _i_ „,I _i_ „,AP 



(29) 



Substituting the corresponding response and correlation functions in Eqs.(20) and (25), we 
obtain the contributions of polar optic phonons as 



7o 



Qp 1 QqC^ 



27r rne* Jq 



{No + 1) (a;^ + no) exp 



k-" 



--[uj^ + noYrtik^) 



+ 



+iVo (cu_L - 1^0 ) exp 



--{uJ^-noYT'Ak±) 



(30) 



and 



dk\k\ 



2V27r m^e* 



(A^o + 1) exp 



dk 



rc{ki_)f{k,) 



A;2 



--{uj + u^ + Qo) r^(k±) 



+ 



+ exp 



-| {u;-u;±- noYr^{k±) 



+ 



+No exp 



--{u; + u;±- Oo) 7-^(^±) 



+ 



+ exp 



^{uj - uj± + QQfT^{kj 



(31) 



the contributions of deformational acoustic phonons as 

1 



A'K\/2'K mfm Jo 



oo p+oo 

i/v-^ I dkMl{k^)p{k,) 



dk\k^ 



(iVfc + 1) {uj± + uk) exp 



~(a;x + M/c)V^2(A;^) 



+ 



+Nk {u!± — uk) exp 



-^{uj±-uk)\^{k^) 



(32) 



and 



1 hD^ 



+00 



87rV27r m^pu Jq 
{Nk + 1) ( exp 



dkzkTc{k±)f{kz) 



■^{uj + u!± + ukf r^{k±) 



+ 



+ exp 



--{(jj -u^-uk) rj/cx) 



+ 



+Nk exp 



--{uj + u;± -uk) r^{k±) 



+ 



+ exp 



{lu -uj± + ukf T^{kA 



and, finally, the contributions of the charged impurities as 



1 eV: 



4„* /■+00 



* ' dk±k^ 



(A;2 + r-Q ^) 



and 



1 



p4 * /'+00 

/ dk^kl 



27rV27r m^eg Jo 



exp 



--{u + u^Yr^ik^) 



+ exp 



(33) 



(34) 



(35) 



In all these formulae we use the notation k = ^/k'j_ + k'^ and uj± = hk\/2m. Nq — 
[exp{—hflo/T) — and Nk — [exp{—huk/T) — 1]'^ are Bose distribution functions for 
the optical and acoustic phonons, respectively. The level of thermal velocity fluctuations, 
{V^), can be determined self- consistently using (V^) = K{ijJc)/2'^q [8]. 



IV. SPIN DYNAMICS 



In the second stage we analyze spin relaxation due to DP interaction between spin and 
electron orbital motion. The corresponding spin interaction Hamiltonian is given by [2,3] 



(36) 



where hb • is the Zeeman splitting Hamiltonian term, /^s is the Bohr magneton, 
and 



with 



(37) 



X. = mV, {V^' - VT) , Hy = rT?Vy {yl - Vl) , ^ mV, {V,' - if) . (38) 
For a quantum well strongly confined in < 001 > direction, Hop can be simplified as [3] 

Hdp = -CTxQxit) - (7yQy{t), (39) 

where 

Qx{t) = XV,{t), Qy{t) = -XVy{t), (40) 

and 

with (p^) = mhiijQ/2 if the confinement is parabolic. 

Considering orbital motion to play the role of an effective heat bath with correlation 
functions given by Eqs.(23) and (24), and employing a second application of the method of 
Refs. [7-9] , we obtain a set of Bloch equations for the average spin projections as 

I M^)) = - K + 5.) (a,(t)) , 
^^{ay{t)) ^{uB + 5y){a^{t))-^^^. 

J,M))-'^^^. (42) 

where Ox^Oy^Uz are the Pauli matrices and a\ = — tanh (fk4;s/2T) is the equilibrium z- 
component of spin, cub — g/J'sB/h, /ib — |e| h/2mQC is the Bohr magneton, and the g^-factor 
is —0.44 for GaAs. It should be emphasized that we determine the relaxation times Ti and 
T2 microscopically as 

Ti = |, T2 = T„ (43) 
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where 

1 2A2 K{lob) 



(44) 



Ts (a;B-a;c)^ + 7o' 

While this result is reminiscent of a formula obtained by Ivchenko for bulk semiconductors 
[10], that differs from the present quantum well case under consideration here in that the 
Ivchenko result has a sum of several Lorentzians and involves numerical constants in place of 
our relaxation rate 70 and numerator function K{uj) which have magnetic field and temper- 
ature dependencies that are explicitly determined by Eq. (29) on the microscopic basis. It is 
evident from Eq. (44), that the effect of magnetic field on spin relaxation is twofold: Firstly, 
there is magnetic field dependence of the fluctuation source correlator, K{ijJb)-i associated 
with electron transitions between energy levels having different spin. Secondly, the difference 
between the frequencies ujb and ujc is involved in the denominator of Eq. (44), which occurs 
because of deviations of electron effective mass and g-factor in semiconductors from the free 
electron values. The temperature dependence can also be separated into two parts, (i) via 
electron thermal fluctuations (this contribution is linear with increasing temperature); and 
(ii) via temperature dependence of the momentum relaxation rate, 70 = 70(7")- 

In the absence of a magnetic field — lob — ^, (Yx) — -f^(0)/27o = T/m), we recover 
Eq. (2) as the zero field limit: 

J_ ^ 4A^ ^ 4A^ J_ ^ 4A^ J_ T ^ a' {plf 1 ^ 

rs B=o 7o 270 ■ ^270^"^/ ^2 7o m 2h^w?eg 70 ' 

wherein our microscopic analysis yields the phenomenological constant Tp of Eqs. (1,2) as 

Tp l/7o- 

Figs. l(a,b) exhibit the dependence of the relaxation rate on the energy of the low- 
est electronic subband of the quantum well (which is given by Eq = huJo/2 for the case of 
harmonic confinement) for various temperatures and magnetic field strengths. It is evident 
from these figures that at low temperature the applied magnetic field suppresses spin re- 
laxation, whereas at higher temperatures the momentum relaxation rate dominates in the 
denominator of Eq. (44) and there is negligible magnetic field dependence in this case. 
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In these calculations we have employed the following set of parameters for a GaAs-based 
quantum well: electron effective mass m — 0.067mo, 1/e* = 0.0069, optical phonon energy 
hflo — 0.035 eV, deformation potential D = 8.6 eV, density p = 5.4 g cm~^, sound velocity 
u — 5 ■ 10^ cm s~^, impurity concentration — 10^^ cm'^, and static permittivity eo = 13. 
It should be emphasized that our microscopic calculations yield quantitative agreement with 
the experimental results of Ref. [6]. The full magnetic field dependence (up to 1000 Gs) 
of the spin relaxation rate is presented in Fig. 2 for fkoo — 0.01 eV. This dependence has 
the well-known [2] Lorentzian shape, evident in Eq. (44), and for a reasonable range of 
magnetic field strength it is not affected by the magnetic field dependence of the velocity 
fluctuations, which are embedded in the function K{u!b)- (However, for the case of dilute 
magnetic semiconductors having a large electron g^-factor, to which our general analysis is 
also applicable, the magnetic field dependence of the function K{u!b) can be of crucial im- 
portance.) The magnetic field dependence for the case of bulk semiconductors (represented 
by a sum of several Lorentzians) was obtained in Ref. [10] and is also reconfirmed by our 
microscopic analysis [11]. 

The temperature dependence of the spin relaxation rate is shown in Fig. 3. The non- 
monotonic behavior is due to the dominance of different scattering mechanisms in different 
temperature ranges. At low temperatures (before the peak), the momentum scattering rate 
is determined by impurities and is almost independent of temperature, and, consequently, 
we have almost linear growth of the spin relaxation rate with increasing temperature, as 
predicted by Eq. (2). However, as temperature further increases, optical phonons become 
the dominant scattering mechanism and, with this, the momentum scattering time becomes 
temperature dependent, resulting in suppression of the relaxation rate. It is important to 
note that in the temperature range in which phonons dominate, the relaxation rate is almost 
independent of the magnetic field, whereas at low temperatures the magnetic field shifts the 
peak and gives rise to deviations from linear behavior. 
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V. CONCLUSIONS 



In conclusion, we have developed a fully microscopic theory of electron-spin relaxation 
in semiconductors by the D'yakonov-Perel' mechanism. We have applied this theory to a 
quantum well structure with a magnetic field in the growth direction. A set of Bloch equa- 
tions for a spin system has been derived with microscopically determined longitudinal, Ti, 
and transverse, T2, relaxation times, which are related as Ti = T2/2 for the quantum well 
structure grown in < 001 > direction [3]. The well-known semiphenomenological expres- 
sion for the spin relaxation rate [3] emerges as the zero-magnetic field limit of our result. 
Furthermore, we have analyzed the dependencies of the electron-spin relaxation rate on the 
energy of the lowest quantum well subband, on magnetic field strength, and on temperature 
obtaining quantitative agreement with the experimental results. 

APPENDIX A 

The derivation of Langevin-like equations for electron velocity operators, Eq.(15) pre- 
sented here follows the approach developed in our previous works [7-9]. We start from 
Heisenberg equations of motion for electron position and velocity operators for the Hamil- 
tonian of Eq.(3), as 

7-1 -J 

jV^it) = [V,{t), H]_ = -u^Vyit) + - E (^kW + t/k) X^{t), 

k 

fVy{t) = ^ [Vy{t), H]_ = u,V,{t) + ^ E ^^2/ (^i^W + ^x) X^t). (Al) 

k 

If the coupling between the electron subsystem and heat bath is weak, or in the case of 
Gaussian statistics of unperturbed heat bath variables, the fully coupled Heisenberg heat 
bath variable is given by [7] 

Qk(t) = Q'^\t) + r citi<^k(t;ti)X_k(ti), (A2) 

J —00 
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where the response function of the phonon heat bath </?k(^; ^i) is defined by Eq.(lO). Substi- 
tuting Eq.(A2) into Eq.(Al), we have to take into account the fact that only the fully coupled 
heat bath variable commutes with the electron variable taken at equal times. Accordingly, 
we perform a symmetrization of both terms of Eq.(A2) with the electron variables, with the 
results 



di 



k ^ 

+ j dhMt;ti)^ [Xk(t),x_k(ti)]+^ , 

jVyit)-u;Mt) -^J2'^y {l [f^k,^kW]+ + ^ [Q^^\t),Xy,{t)\^ + 



* 1 

dtiipi,{t;ti)-[Xk{t),X^k{ti)] 

oo ^ 



+ 



(A3) 



To eliminate the unperturbed phonon/impurity variables we employ the Furutsu-Novikov 
theorem [12]: 



Q'^\t),x^{t) 



dtiM]^{t] ti 



(A4) 



where the correlation function Mk(t;ti) of free phonon variables is given by Eq.(ll) and 
S/SQ^'^^ti) is the functional derivative with respect to the uncoupled heat bath variable 
Q^^\ti). Eq.(A4) can be derived by an application of the Wick theorem, assuming the 
operator Xk(t) to be a functional of {<5*^°-*(ti)} with ti < t (See Ref. [7]). This expression 
is exact only for Gaussian statistics of the variables Q^^\t), and it also can be applied in 
the case of weak couphng. For the case of strong couphng, Eq.(A4) requires modification to 
include functional derivatives of all orders. The functional derivative on the right-hand side 
of Eq. (A4) is proportional to the commutator [7] in the form 

5X^{t) 



with the following result 



h 



[^k(t),^-k(ti)]_)e(t-ti). 



(A5) 



Q^^\t),x^{t) 



dtiMk(t;ti)( -[Xk(t),X_k(ti)]. 



(A6) 
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An analogous treatment of the electron-impurity correlation term yields 

[C/k(i),Xk(i)]+^ = dti$k [Xk(i),X_k(ti)]_^ , (A7) 

where $k is the impurity potential correlation function given by Eq.(13). We also introduce 
fluctuation source operators defined as 

k ^ 

- dii(Mk(i;ii) + $k)^[^kW,^-k(ii)]_^ (A8) 
with zero average, {^x,y{t)) — 0, and collision terms defined as 

G,,y[yM.Vy{t)] ^--Y^i^y f c?iif(Mk(i;ii) + $k)^[^kW,^-k(^i)]_ + 

^ k J -oo \ 

+V^.{t■M)\ [Xk(i),X_k(ii)]+) , (A9) 
resulting in the Langevin-like equations of Eq.(15). 
APPENDIX B 

In this appendix we determine the commutators [Xk(t), X_k(^i)]-|- involved in Eqs.(16), 
(A8), and (A9). We assume that there is finite correlation time, Tc, of the electron-phonon 
interaction (which will be obtained self- consistently). For the case of weak coupling, the 
commutators can be calculated using the free, uncoupled evolution of the electron velocity 
operators during Tc , as given by 

At) = + 14(ti) " '''^ - ^ " ""^'"-'^ " (Bl) 
y[t) = y{ti) + Vx{ti) h Vy[ti) . 

This allows to us determine the following commutator 

[^k±r^(^), -^kxrx(^l)]_ = sin(cUc(t - ti)), (B2) 

where = |k^| = y^k^~+~k^ is the magnitude of the transverse wave vector. With the 
operator equation (Baker-Campbell-Hausdorff identity) 
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(which is vahd, when A, is a c-number) we obtain 

7 9 f flip' \ 

[Xk(t),X_k(ti)]_ = ^/2(A;,)exp{zkx(r^(t) -r^(ti))}sin sin(u;,(t - t^) j , (B4) 

- [Xk(i),X_k(ii)]+ = exp{ik^ {v^{t) - r^(ii))}cos ( ^-^ sin(a;,(t - t^))\ . 

We assume the coordinate fluctuations to be approximately Gaussian (see Ref. [8] for the 
corresponding discussion) and, consequently, obtain 

(cxp (r^(t) - r^(ti))}) = exp |-i ((k^ (r^(t) - r^(ti)))^>} (B5) 

and 

exp{ik^ {v^{t) - r^(ti))} - (exp{ik^ {r^{t) - r^(ti))}) (B6) 
~ ik^ (r^(t) - r^(ti)) (exp{zk^ (r^(t) - r^(ti))}) . 

With these simplifications we have 

G.AyS):Vy{t)] - (Gx,.[K(t);v;(t)]) = (B7) 

= ih,yfiK) J dh I (Mk(i; ii) + $k) ^ sin (^^-^ sin(a;,(i - ii)) j + 

+ipk{t; ti) cos (^^^ sin(u;e(^ - ^i))^ | {r±{t) - r^(ti)) (exp {ik^_ {r^{t) - r^(ii))}) . 

We make the further assumption that the correlation time of the electron-phonon interaction, 
Tc, is much less than the period of the cyclotron oscillations, i.e. lOcTc « 1 (the same 
approximation was used in Ref. [10]). This is reasonable for semiconductors at moderate 
magnetic fields and not too low temperatures, and it leads to 

k± {r±{t) - r^{ti)) ^ k^V^{t){t - h) + kyVy{t){t - h), (B8) 

-sm(ujJt-ti)) ^ -(t-ti), (B9) 

16 



and 

\2 



(exp{ik^ {v^{t) - r^im) = exp |-|^| , (BIO) 
where 

re-'(k±) = kl (V^it)) + kl {V^{t)) + k^ky Vy{t)]^) . (Bll) 

The resulting simphfied stochastic equations for the fluctuating velocity components take 
the form given by Eq.(19). 
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FIGURES 

FIG. 1. Dependencies of the relaxation rate on the energy of the lowest level of the quantum 

well; (a) for temperature T = 40K and two magnetic field strengths, (b) for temperature T = 300K 
(the same curve for both magnetic field strengths). 

FIG. 2. Dependence of the relaxation rate on magnetic field strength for various temperatures. 

FIG. 3. Temperature dependence of the relaxation rate for various magnetic field strengths. 
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